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US Report to the UNFCCC 

Annual Inventory of Greenhouse Gas Emissions and Sinks 

In 2010: 

US emissions from fossil fuels: 1,556 Tg C 

US forest sink: 251 +/- 80 Tg C or 16% of fossil fuel emissions 

Forest sink estimate from FIA (periodic sampling of >150,000 plots) Sold Standard" 
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Reasons why a forests would be carbon sinks 
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C0 2 , Precipitation, Solar Radiation 

Nitrogen, Temperature 


NEP = NPP - RH ( = GPP - RE) 

***Regrowth sink: reduced carbon pools 
Enhancement: increasing NPP 

Other processes less understood: decreases in mortality, respiration attenuation 


Forest Sector Carbon Exchanges 



Translate A C stocks intoTVEEto diagnose regrowth vs enhancement 

NEP stock change = AC s,ocks + Harvest+ Fire Emissions 

Includes all effects including age and enhancement 

NEP 'regrowth = NPP-RH 

Model accounts for age only (same as Canadian approach for UNFCCC Reporting) 


Our Approach: 

-Calibrate CASA with FIA Forest Age-Biomass data specific to forest type, productivity level and region 


-Develop a spatially explicit map of forest age for CONUS (from FIA and Landsat based VCT, Goward et al) 


-Assign carbon fluxes based on age of forest 
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Flux Trajectories 



Unique trajectories for regional forest types, productivity classes, and climate settings 
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Uncertainty formally analyzed with Monte Carlo fitting 

biomass sampling error from FIA (+/- 10 to 100%) 
volume to carbon conversion (+/- 7%) 


Williams et al. (2012) Glob Biogeochem Cyc 




Flux Trajectories Within Regions 
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Williams et al. (2012) Glob Biogeochem Cyc 
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Age Structures Within Regions 
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Regrowth Uptake 
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Williams et al. (2012) Glob Biogeochem Cyc 





Regrowth Uptake Uncertainty 



c u 


NEP regrowth modeling 
Standard Deviation 
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NEP from Stock Changes 
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stock change 

= AForest Stock + Harvest Flux + Fire Flux 
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Why are they so different? 

^ E^stock change ^^^regrowth modeling 



310-160= 150 Tg C a 1 

Stock Change NEP is 2 X larger than Regrowth sink 
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What does the Landsat VCT disturbance history (Goward et al.) 
tell us about forest age and regrowth? 
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■ generate detailed maps of forest NEP, Biomass, etc. sources 
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Age Structure of Forests Estimated Using 
Landsat VCT 1985-2005, FIA < 1985 

Introducing fine scale disturbance data reduced the estimated CONUS NEP 

Higher disturbance rates especially in regions with large forest industry produce 
younger forests with lower NEP 


*Makes the discrepancy between Stock Change NEP and Regrowth NEP larger 
Stock Change NEP 2.6 X larger than Regrowth NEP 


East 


AREA NEP RS SE NEP FIA SE NPP RS NPP FIA Wood R 

(le9m 2 ) (TgC/yr) (TgC/yr) (TgC/yr) (TgC/yr) (TgC/yr) (TgC/yr) (TgC) 


212 

0 

355 

420 

192 

493 

127 

202 


10 


1.29 


12 


1.29 


85 88 


0 

0 

0 

0 

0 

0 

22 

3.47 

30 

3.47 

308 

317 

28 

4.27 

40 

4.27 

316 

347 

4 

1.81 

7 

1.81 

65 

73 

10 

5.51 

11 

5.51 

87 

88 

8 

2.81 

13 

2.81 

89 

97 

10 

2.93 

18 

2.93 

167 

171 

120 

27.88 

164 

27.88 

1377 

1451 


706 

0 


851 

0 


834 1236 


FIA <25yr RS <25yr FIA 

<5yr RS 

<5yr FIA 

5 

1.46721854 

3 

:) % 

% 

% 

% 



23 

16 

4 

3 


0.95 


0 

0 

0 

0 


0 


51 

39 

12 

8 


2.27 


45 

37 

10 

8 


1.96 West 

34 

21 

8 

5 


1.55 

1.416992 

11 

26 

1 

7 

0 


1.2 


28 

11 

9 

2 

1.29 


45 

19 

8 

4 


1.9 


30 

17 

7 

4 


1.34 



Total 


2339 


11094 


13015 







Sensitivity Analyses 


Errors in the representation of biomass accumulation with stand age could cause 
large errors in estimated fluxes 

1) wood turnover has a large influence on rate of biomass increase and is poorly 
constrained 


2) stand thinning (selective logging) may lower biomass in the inventory without 


resetting age 

3) age reporting may be biased 

Original Results 

Wood Turnover Increased 10% 

Biomass Increased 10% 

Stand Ace -5 years 

Stand Age +5 years 
w * 
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Biomass would have to be 70% higher to produce the NEP stockchange 


Williams et al. (2012) Glob Biogeochem Cyc 


McKinley et al 'll from Woodbury 
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Recent Studies presenting Evidence of Growth Enhancement 


Methods Employed: Inventory, Modeling, Eddy Covariance, Tree Ring 


Approaches 
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Plots Show Accelerated Growth 



Present-day growth rates exceed those expected from chronosequences of 

biomass regrowth with forest age 
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McMahon et al (2010) PNAS 
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Analysis of individual FIA Plot data for Eastern US Forests 

(>1000 Plots in Each State) 

Construct age dependent biomass trajectories for major forest types (Regrowth) 

Estimate the change in biomass of individual plots between the 2 most recent 
inventories (Stock Change) 



FIA Biomass vs. age 




Biomass change rate 
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Median, 25th %tile, & 75th %tile values of 


biomass change from plot remeasurement 
(~ 5 year interval), in 10-year cohorts 
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Derivative of Biomass vs age curve 
(dashed = 75th %tile, solid = median) 
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85 year old forests 
^ Stock Change= 3 MgC/ha/yr (300 gC/m 2 /yr) 


Enhancement = 3 ° — = 2 
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State - Forest type 


Age-weighted Growth Enhancement by FIA Type/Region 
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Evidence for Enhancement in Eddy Covariance Measurements 

25 forest sites US and Southern Canada. 
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Most sites show positive trends in GPP 
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Data courtesy of K. Schaefer 




Inversions can't yet reconcile magnitude and distribution of Forest Sinks 
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(a) NEP from Age-Accumulation Modeling 


Mean of Inversions 


(b) NEP from Stock Change 






0.4% yr' 1 increasing Woody NBP, US 


Growth Enhancement Scenario: 

Increase Wood NPP at 0.4% / yr for 100 years 
Produces a 90 Tg C / yr Sink 



0.4% increasing Woody NPP, US 





Conclusions 


Modeling disturbance and regrowth effects on forest NEP accounts for 
only half of the carbon sink seen in US stock changes 

The additional sink is likely due to growth enhancements (climate 
trends, C02 fertilization, nitrogen fertilization, irrigation) 

Enhancement is concentrated in regions with active management and 
younger and more productive forests (SE, SC, PNW, PSW) 

Enhancement is supported by a number of lines of evidence. 

The mechanisms driving the enhancement are unclear. 


Note the Size of the enhancement exceeds expected C0 2 fertilization effects so positive 
feedbacks with other factors are required 


Current and Future Work 



♦ Wall-to-Wall CONUS VCT Products (PI: Goward) 

♦ Expand spatial domain to Canada and Alaska 

♦ Improved representation of disturbance type (fire, harvest, insect driven mortality) 

♦ Derive finer spatial scale calibrations from FIA (state?) and account for thinning. 

♦ Address the issue of what if anything is causing large scale growth enhancement 

♦ Evaluate ability to detect/quantify forest source/sinks from observed atmospheric 

C0 2 variability (Kawa et al.) 




Detecting and Quantifying Vegetation Responses to Climate Variability Using NDVI 


Jim Collatz, Fanwei Zeng, Alvaro Ivanoff, Jorge Pinzon 

In support of and supported by several projects (GFED, CMS, Kawa) 

Rational: NDVI is used to drive the carbon cycle in diagnostic models like CASA. It is assumed 
that some part of NDVI variability represents real responses of vegetation to climate 
in these models. To do a better job at predicting interannual variability in the carbon cycle we 
need to understand the sources of variability in the models. 

How much of the NDVI variability is caused by real climate response of the vegetation? 

Data Sets: 

NDVI: GIMMSg AVHRR, MODIS Terra, MODIS Aqua 
Temperature: GISS and CRU 
Precipitation: GPCP, TRMM 

Other: 

MODIS Snow Cover, MODIS Aerosol OD, MODIS Cloud Cover, GISS Solar Radiation 


Weakness in AVHRR addressed by MODIS: degradation of calibration, drift, other BRDF, narrower bands (corrections for 
water vapor and aerosols), resolution. 

Advantages of AVHRR (esp GIMMSg): long timer period (more robust detection of climate signal), long experience with data 
(> 25 years) 


GIMMSg NDVI captures same climate (and noise) driven variability as MODIS 



Significant correlations (p < 0.05) between GIMMS and 
(a) MODIS Aqua monthly NDVI anomalies, 2003-2010 



Significant correlations (p < 0.05) between MODIS Terra and 
O 3 ) Aqua monthly NDVI anomalies 2003-2010 




Significant correlations (p < 0.05) between GIMMS monthly Significant correlations (p < 0.05) between MODIS monthly NDVI 

NDVI anomalies and antecedent TRMM precipitation, anomalies and antecedent TRMM precipitation, 2000-2010 





Figure 4. Global maps of month-specific correlation coefficients between GIMMS NDVI and cumulative GPCP precipitation 
(1-6 month lead) for 1982-2010 (n = 29). Only correlations significant at 95% confidence level are shown. (Spatial 
resolution: l°xl°. NH: Northern Hemisphere; SH: Southern Hemisphere) 




r: 1962*2010 Jan GMMS NDVI - GPCP cum precipitation (lead time = 1 to 6 months. p«0 0$) r: 1982-2010 Feb G*tMS NDVI - GPCP cum precipitation (lead time s 1 to 6 months, p<0 05) 


r: 1962*2010 Mar GIMMS NOVI - GPCP cum precipitation Head time = 1 to 6 months. p<0 05) r 1962*2010 Apr GIMMS NOVI • GPCP cum precipitation (lead time = 1 to 6 months. p-«0 05) 
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r 1962*2010 May GJMMS NDVI * GPCP cum precipitation (lead time s 1 to 6 months, p<£ 05) r 1982*2010 Jun GIMMS NOVI - GPCP cum precipitation (lead time = 1 to 6 months. p«0 05) 
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Figure 4. The correlation coefficients between GIMMS NDVI and cumulative GPCP precipitation (1-6 month lead) the 
month of July, 1982-2010 (n = 29). Only correlations significant at 95% confidence level are shown. 

r: 1962-2010 Jul G4AMS NDVI - GPCP cum precipitation (lead time = 1 to 6 months. p<0 OS) 



pixel % that show significant positive correlation (95% confidence level) and average correlation coefficient aggregated into 
certain vegetation classes and regions. Regions and classes are examples where high correlations occurred. 




Correlations show that NDVI captures the seasonal dynamics of ecohydrology 


Figure 5. Regional maps of month-specific GIMMS NDVI and GPCP precipitation correlations with varying lags (n = 29). 
(a) North America: 


Lag - I month Lag 2 months Lag - 3 months Lag 4 months Lag $ months lag 6 months 



Early in the North American growing season vegetation growth is sensitive to 
cumulative antecedent precipitation from winter 

Late in the growing season vegetation is more sensitive to recent precipitation, 
similar behavior in Africa and Australia. 


Significant (p < 0.05) correlations between temperature and NDVI anomalies 


Figure 7. Global maps oltnotuh- specific correlation coefficients between GIMMS NDVI and GISS temperature anomalies with lag 
0 month for 1982-2010 (n " 29), after removal of pixel* with significant negative correlation between GIMMS NDVI and MODIS 
Terra snow cover Only correlations significant at 95** confidence level arc shown. (Spatial resolution: 0.5' *0.5*. NH Northern 
Hemisphere: SH Southern Hemisphere) 
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Note: negative correlations are associated with negative correlations between precipitation and temperature 
-water limited conditions 


NDVI and temperature anomalies are correlated with the onset of spring in northern latitudes 


The variance in NDVI explained by "real" vegetation responses to climate 
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sum <82*10. Jul>: positive CIMMS NOVI • cum GPCP and GISS <p<0 05, no snow) r 


sum <82*10, Sep): positive GIMMS NOVI - cum GPCP and GISS <p«0O5. no snow 




sum (82*10, Aug): positive GIMMS NOVI • cum GPCP and GISS <p<0 05. no snow) 


sum <82*10, Oct): positive GIMMS NOVI • cum GPCP and GISS <p*0 05. no snow 



sum <82*10, Nov) positive GIMMS NDVI * cum GPCP and GISS (p<0 05. no snow) 
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Variance in NDVI caused by noise 


Variance in NDVI explained by current month precip, snow, aerosols, and clouds 

Example: September 

**Note all correlations are negative 

R 2 (00-10. Scpi GIMMS - GPCP, snow, aerosol and cloud (load = 0 mon lor all. p«0 05) 
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NDVI interannual monthly variance explained by 


"real" climate responses 

— — - Interference from snow, precipitation, clouds and aerosols 
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Atmospheric C0 2 inversion 
June-August 2003, Peters etal.,GCB'10 
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